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ABSTRACT 



in 

Context. The more massive counterparts of T Tauri stars, Herbig Ae/Be stars, are known to vary in a complex way with no variability 
i-C ■ mechanism clearly identified. 

Aims. We attempt to characterize the optical variability of HD 37806 (MWC 120) on time scales ranging between minutes and several 
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years. 

Methods. A continuous, one-minute resolution, 21 day-long sequence of MOST (Microvariability & Oscillations of STars) satellite 
observations has been analyzed using wavelet, scalegram and dispersion analysis tools. The MOST data have been augmented by 
sparse observations over 9 seasons from ASAS (All Sky Automated Survey), by previously non-analyzed ESO (European Southern 
Observatory) data partly covering 3 seasons and by archival measurements dating back half a century ago. 

Results. Mutually superimposed flares or accretion instabilities grow in size from about 0.0003 of the mean flux on a time scale 
of minutes to a peak-to-peak range of < 0.05 on a time scale of a few years. The resulting variability has properties of stochastic 
"red" noise, whose self-similar characteristics are very similar to those observed in cataclysmic binary stars, but with much longer 
^\ ' characteristic time scales of hours to days (rather than minutes) and with amplitudes which appear to cease growing in size on time 

^\ , scales of tens of years. In addition to chaotic brightness variations combined with stochastic noise, the MOST data show a weakly 

I/") . denned cyclic signal with a period of about 1.5 days, which may correspond to the rotation of the star. 
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1. Introduction as A2Vpe or B9 Ve+sh (Swings & Strove [19431 Guetter [19811 

The et al. |1994| l although sometimes it is still quoted as AO. 

HD 37806 (SAO 132452, RA(2000) = 05:41:02.3, DEC(2000) HD 37806 was one of the first << ra dio" sta rs identified (Feldman 

= -02:43:01, V - 7.95) is an early star in the field of Sigma Ori. et a j |1973[ ) 
It is frequently referred to as MWC 120, the designation in the 

b Merrill & Burwell Catalog ([I933T I of early type stars with emis- The age and cluster membership of HD 37806 are not well 

■ - ■ sio n lines , which it received following the discovery of Merrill et known, although the star does appear very young. While it is 

al. (fT932l . HD 37806 is a Herbig Ae/Be-type object, i.e. a young only 35 arcmin from Sigma Orionis, it was not included among 

star which is more massive (typically 1 .5 - 5 M ) than solar and *e members of the young cluster carry ing the name of this star 

sub-solar mass T Tauri stars. The review of Waters & Waelkens «i the study by Sherry et al. ( 2004 , 2008 ). The distance estimated 

( 1998) summarizes the properties of Herbig Ae/Be-type stars; a for the cluster, 420 + 30 pc, is consistent with the Hipparcos 

ca talogu e of Herbig Ae/Be-type objects was published by The et new-reduction parallax 2.05 + 0.9 1 mas (van Leeuwen |2007|l. 

al. rn?94T >. The MK spectral type of HD 37806 is variously given The star almost certainly belongs to the Orion OB lb association 

(of which the Sigma Ori cluster may be a part), whose age is 

I „ j Z . c ~T~ i.r^c-r . j- c a estimated at 2 - 5 Myr (Brown et al. 1994) and for which the 

Based on data from the MOST satellite, a Canadian Space Agency TT . J . ,.„,i9n — ; — i , i, „» A , 

mission jointly operated by Dynacon Inc., the University of Toronto Hipparcos group distance is 350!™ pc In Brown et al. 0994), 

Institute for Aerospace Studies and the University of British Columbia, the distance to HD 37806 was estimated at about 300 pc, assum- 

with the assistance of the University of Vienna, and on data from the All "lg V = ISA, My = +0.2, and a moder ately large Ay = 0.35; 

Sky Automated Survey (ASAS) conducted by the Warsaw University for the Sigma Ori cluster, Sherry et al. (120041 [2008 ) estimate 

Observatory, Warsaw, Poland at the Las Campanas Observatory, Chile. Ay at about 0.1 to 0.2. The distance, the possible membership of 
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HD 37806 to the Sigma Ori cluster and the amount of interstellar 
absorption all appear to be very uncertain at this point. 

The main thrust of this paper is limited to utilizing the con- 
tinuous, 21 day long observations by the MOST satellite for 
characterization of the optical variability of HD 37806. As aux- 
iliary data, we use lower accuracy ESO and ASAS observations, 
which extend the time coverage to several years, and archival ob- 
servations extending the time baseline to half a century. The data 
analysis tools are similar to those used for continuous MOST 
satellite observations of the T Tauri type star TW Hya (Rucinski 
et al. 120081 . 

The literature on diversified studies of HD 37806 in various 
spectral regions is extensive, where the star is usually identified 
by its more popular alias MWC 120. Among the recent studies, 
Mottram et al. d20071 > pointed out that Herbig Ae stars (A0-A2), 
such as HD 37806, appear to be more similar in terms of the ac- 
cretion processes producing the emission spectra to T Tauri stars 
than to earlier (B0-B3) Herbig Be stars. The star has also been 
the subject of extensive spectroscopic, polarimetric and spectro- 
polarimetric investigations (Vink et al. 120051 1. Of relevance to 
variability studies of HD 37806 is the finding by Harrington 
& Kuhn (2009) that the Ha line shows strong morphological 
changes on a time scale of 3 years. 

Very recently, using a novel method of spectro-astrometric 
signatures, Wheelwright et al. (120101 1 found that HD 37806 is a 
binary system, similarly to the majority of surveyed 45 Herbig 
Ae/Be stars. This is a very important discovery which may pro- 
vide the key explanation of the variability characteristics which 
we see in our MOST data; we return to this subject in Sections|4] 
and |5] 



2. Optical variability: Previous photometric results 

2.1. Archival data 

HD 37806 was observed in the Johnson V and Stromgren y 
(tra nsform ed to V) bands in 1959 - 1962 and 1968 by Hardie et 
al. d 196411 an d Warren & Hesser (11 9771 1. Additionally, de Winter 
et al. J200TT > observed the star in 1986 and 1992. Although all 
these were few (2 and 5) scattered observations, they give an 
interesting constraint on the long term variability of the star 
(Section©. 

2.2. Hipparcos observations 

Hipparcos satellite observations of 44 Herbig Ae/Be stars, ob- 
tained over 37 months of the mission in the Hp band, were 
analyzed by van den Ancker (1998 ). Typically a hundred time- 
dispersed Hp magnitude determinations per star were obtained 
with accuracy of 0.008 - 0.015 mag. In this study, HD 37806 
was observed to vary by about 0.07 mag. However, this study 
was obviously unable to define detailed variability characteris- 
tics of the star. 

We looked at the Hipparcos data again. The data consist of 
135 Hp observations and 209 Vp and Bj observations of the 
Tycho project, the latter with rather large uncertainties of typ- 
ically 0.06 - 0.12 mag. All observations tend to be strongly 
bunched in time as a result of the sky scanning pattern of the 
satellite and appear in only a few groups over the mission dura- 
tion. The strong bunching would not be a detrimental obstacle 
for use in the long time-scale variability considerations, as in 
Section |U if the main-project data were in the V rather than Hp 
system. 
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Fig. 1. The ESO light curve of HD 37806 in 3 seasons of 1985 - 
1986, expressed in units of the mean flux. The central grouping 
of points corresponds to one Chilean observing season; the two 
other seasons were poorly covered. The HJD offset used in this 
figure is the same as for the remaining plots in this paper, hence 
the negative values of time. 



In the end, we used only the Tycho-2 (H0g et al. 2000) 
well-calibrated, single data point corresponding to the average 
V = 7.934 + 0.016 (note that there exists no Tycho-2 epoch pho- 
tometry). It has been added to the remaining archival observa- 
tions discussed later on (Section |4| in an attempt to relate vari- 
ability seen in by us to its possible long-term continuation. 

2.3. ESO observations 

Four-color Stromgren uvby observations of HD 37806 were con- 
ducted at the European Southern Observatory (ESO) in 1985 - 
1986 by several observers (Manfroid et al. 1991); the observa- 
tions do not appear to have been analyzed so far. Altogether 55 
observations in the y (= V) band were obtained, with 12 dur- 
ing the first season of January - February 1985, 41 observations 
during the November 1985 - March 1986 season and with 2 ob- 
servations in September 1986. The typical interval between ob- 
servations was 2 days and there were no multiple observations 
taken on the same night. The observations were done differen- 
tially and were relatively accurate with typical measurement er- 
rors of 0.007 mag. 

The ESO observations of HD 37806 are shown in Figure Q] 
The observations do not indicate any regularity in its variability, 
but rather random noise, perhaps with increasing amplitudes for 
time scales ranging from single days to hundreds of days. We 
discuss the optical flux scatter in Section [3~3l where we inter- 
pret the variations in terms of stochastic noise characterized by 
c\ = cr 2 , - <j 2 R , where the observational (instrumental) noise is 
described by the contribution <x R = 0.007. 

2.4. ASAS observations 

ASAS, the All Sk y Automated Survey (Poimanski [19971 [20021 
2004, 120051 12006^1). is a long term project dedicated to detection 
and monitoring of variability of bright stars using small tele- 
scopes. It has been run for a decade by Warsaw University at 
Las Campanas Observatory in Chile using 7 cm telescopes pro- 
viding photometry in the 8-13 magnitude range. About three 

1 See: http://www.astrouw.edu.pl/asas/ 
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Fig. 3. The MOST light curve of HD 37806, covering 21 days, expressed in units of the mean flux, with one-minute sampling. The 
first half of the run experienced gaps of about 30% of each 103 minute satellite orbit due to a high background level. The data set has 
been used at the full one-minute resolution as well as at the uniform, MOST-orbit time sampling necessary for the wavelet analysis 
(see Section [3~!2T ) 
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Fig. 2. The AS AS light curve of HD 37806 expressed in units of 
the mean flux. The groupings of points correspond to individual 
observing seasons. Note the different horizontal scale compared 
with Figure Q] 



ing the large field of view, extinction anomalies - random er- 
rors when individual frames are compared, as in variability stud- 
ies. Also, the inexpensive (engineering grade CCD) detector had 
many bad and warm pixels which were observed to vary their re- 
sponse in time and were sometimes not immediately recognized 
in automatic reductions. These issues, coupled with the rather 
small variability amplitude for HD 37806 forced re-reductions 
of the data for this star which are slightly different from those 
available from the automatic on-line ASAS data system. 

The observed brightness variations of HD 37806 are dom- 
inated by random stellar flux fluctuations, enlarged by the ob- 
servational noise. In the analysis of the brightness fluctuations 
in Section 13.31 we assume observational random noise of <tr 
().() 1 7. The two well defined brightness decreases (at time 3050 
and 4100 in Figure |2| appear genuine as we could not identify 
any instrumental effects to produce them; they may be mani- 
festations of sudden obscuration brightness drops observed for 
the UX Ori-type subclass of Herbig Ae/Be stars (Waters & 
Waelkens 1998). 



fourths of the sky (-90° < 6 < +28°) is covered on a regular 
basis. 

The ASAS observations of HD 37806, 424 in number, were 
obtained in the V filter over 9 seasons, from February 2001 to 
December 2009. They are shown in Figure [2] The typical ran- 
dom errors are somewhat difficult to characterize and were typi- 
cally at the level of 0.015 - 0.02 mag. The errors include a major 
component resulting from reductions of individual frames to the 
V band, which are done relative to all available Tycho-2 cata- 
logue stars within the field of view of the telescope. While this 
approach assures the photometrically differential nature of the 
data (with possible systematic errors within the frame), it may 
generate - due to environmental factors such as clouds cross- 



3. MOST satellite observations 

3.1. The data 

Originally planned for two years, the MOST satellite has been 
observing variability of stars highly successfully since 2003. The 
optical system of the satellite consists of a Rumak-Maksutov f/6 
15 cm reflecting telescope. The custom broad-band filter covers 
the spectral range of 380 - 700 nm with effective wavelength 
located close to Johnson's V band. The pre-launch characteris- 
tics of the mission are described by Walker et al. (2003) and the 
initial post-launch performance by Matthews et al. ([2004). 

HD 37806 was observed by MOST in Guide Star mode for 
21 days, during observations of Sigma Orionis and its surround- 
ing stars in November 2007 (Townsend et al.; in preparation). 
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The data were reduced with other stars of this campaign using 
the MOST Guide Star Photometry reduction method (Hareter et 
al. 12008 ) and expressed as added counts per one minute inter- 
val. In this paper we utilize only the flux values corresponding 
to counts normalized to their mean value. The mean count rates 
are not relevant for error estimation because of the highly non- 
Poissonian noise characteristics resulting from addition of fre- 
quent CCD readouts (typically every second) of the Guide Star 
mode observations. 

Continuous, uninterrupted, one-minute interval observations 
were done during the second half of the run. Thus, a full analy- 
sis at this temporal resolution could be done over a shorter inter- 
val of 9.5 days only (Section [3.3b . During the first half of the 
run, the observations had stray-light elimination gaps of typ- 
ically 29 minutes per 103 minute MOST orbit. In this paper, 
we used both the whole dataset sampled at one minute as well 
as "normal points" obtained by median averaging of naturally 
formed groups of observations at MOST-orbit intervals of 0.07 
d. The orbital normal points were used for the wavelet analysis 
( Section [JSt which requires data sampled at strictly equidistant 
points in time. 

Well-defined fluctuations of a few percent from the mean 
flux are visible in the MOST light curve (Figure [3]). The fluc- 
tuations appear to be reminiscent of flares or irregular brighten- 
ing events observed in cataclysmic and nova-like variables (CV) 
and which are explained as stochastic variability resulting from 
superposition of overlapping accretion instability events in the 
disk surrounding the star. 

Fourier analysis of the MOST data was done for all avail- 
able data sampled at one minute intervals by fitting cosine-sine 
pairs with frequencies incremented in small steps to cover the 
whole accessible frequency range. This corresponded to a range 
extending from the lowest accessible 0.05 cycles per day (c/d) to 
several hundred cycles per day. The amplitude spectrum, shown 
in Figure|4]in logarithmic units, is dominated by low-frequency 
variability at frequencies below a few cycles per day. Usually, 
due to stray-light effects, passages through the South Atlantic 
Anomaly, etc. the MOST orbit frequency of 14.2 c/d, along with 
a frequency of 1 c/d appear in such spectra but both are absent 
in the reduced material, attesting to good processing of the raw 
data. The low frequency content has led us to use in the wavelet 
analysis (Section [3.2b the data averaged into points sampled at 
the MOST-orbit frequency, i.e. below 7 c/d. 

In addition to the Fourier component amplitudes, Figure [4] 
gives the mean standard errors for the amplitudes, estimated us- 
ing a bootstrap experiment of multiply re-sampled data. The er- 
rors start to become significant for frequencies higher than a few 
c/d; their value was used to estimate the average mean random 
error for the MOST data at cr R = 0.0028 of the mean flux. 

No coherent, periodic signal is visible in the amplitude spec- 
trum except for a somewhat complex feature at about 0.65 c/d. 
The overall spectrum appears to be that of red stochastic noise, 
possibly even "redder" (with stronger low frequencies) than the 
common "flicker noise" (amplitudes scaling as a oc 1/ Jf; Press 
119781 ); the slope of the amplitude spectrum may in fact show a 
change around the discrete spectral feature at =* 0.65 c/d or the 
period about 1.5 days (see Figure|4]i. 

3.2. Wavelet transforms 

The data for the wavelet transform must be arranged into a 
strictly equidistant time sequence to accommodate progressive 
extension of the temporal scale by integer multiples of the ba- 
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Fig. 4. The amplitude spectrum of the MOST variability data 
for HD 37806 is shown by the continuous line. The dotted line 
gives the Fourier amplitude errors estimated by a bootstrap re- 
sampling experiment. The broken lines, shown to guide the eye, 
give the slopes of a oc 1/f and a oc 1 / Jf ; the latter corresponds 
to "flicker noise" (Press [T978l . 
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Fig. 5. The Morlet-6 wavelet (sine-cosine pair extending over 6 
cycles and confined by a Gaussian) applied to the HD 37806 data 
sampled at the MOST orbital period of 0.07046 d. The light gray 
areas correspond to stronger power. Edge effects influence the 
wavelet transform beyond the broken lines. Note that the phase 
information for this complex-function transform has been disre- 
garded by using the power of the transform. The missing orbital 
data were interpolated at two points marked by tick marks along 
the lower horizontal axis, but no detrimental effects appear to be 
present there. 



sic sampling interval. For MOST data, the most obvious such 
sampling interval is the duration of one satellite orbit around the 
Earth. Indeed, the majority of instrumental and stray-light prob- 
lems usually manifest themselves with a period of 0.07046 d. 
The median flux values for HD 37806 from individual satellite 
orbits were obtained from typically 70 to 100 data points per me- 
dian in the first half of the run (when orbital gaps were present) 
and up to 102 points (minus outliers and glitches) per median in 
the second half (the last 9.5 days) of the run. The resulting data 
were additionally spline interpolated - with very slight shifts of 
less than a minute - into a strictly equidistant time grid with steps 
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Fig. 6. The same as in Figure [3] but for the DOG-2 wavelet (the 
negative second derivative of a Gaussian). Note that the shortest 
coherent noise structures seem to appear at time scales of about 
0.25 d. 



of 0.07046 d. Two gaps in the mean data were present on HJD - 
2,450,000 days 4427 and 4437; they were filled by spline inter- 
polated points to obtain a uniform sequence of 301 data points. 

Wavelet analysis is normally done with time-scale incre- 
mented as a power of 2, i.e. in 2, 4, 8, 16... data sampling inter- 
vals (here 0.07046 d each). The same sequence was used here, 
but - for convenience of the subsequent transform gray-scale dis- 
plays - the scales were later interpolated into a linear sequence 
extending from 0.14 to 4 days; for still longer time scales, edge 
effects start becoming important for the MOST data. 

Various types of wavelet functions are in widespread use 
in data processing (Torrence & Compo 1998). Selection of a 
particular type depends on the type of variability. The simplest 
and most often used are the Morlet wavelet of a sine-cosine 
pair confined in a Gaussian wave "packet", which is useful for 
detection of short-lasting or evolving periodic events, and the 
Derivative of a Gaussian (DOG) wavelet, in its lowest-order ver- 
sion (the second derivative, DOG-2), which is sometimes called 
the Mexican Hat wavelet and which is useful for characterization 
of chaotic, localized events. We used both because variability of 
HD 37806 appears to be dominated by random, chaotically oc- 
curring flares, yet the lesson of TW Hya (Rucinski et al. 2008 ) 
taught us that such stars may show hidden, evolving periodici- 
ties. For the Morlet wavelet, we used the 6-th order version be- 
cause it was found in the analysis of TW Hya that the more com- 
monly used Morlet-5 did not reproduce the Fourier frequencies 
sufficiently well (and this has been confirmed again here); qual- 
itatively, the results for both Morlet transforms are very similar. 

The Morlet-6 transform (Figure [5]) confirms the presence of 
a quasi-periodic signal with a period of about 1 .5 d, which was 
already noted in the Fourier analysis (Figure @}. In the wavelet 
transform, the 1.5 d signature appeared twice during the three 
weeks of the observations, with wavelet amplitudes 0.027 and 
0.040 of the mean flux, which can be compared with the am- 
plitude of a bit less than 0.003 in the whole, "diluted-in-time" 
Fourier spectrum. In Section|4]we interpret the 1.5 day period as 
related to the rotation of the star with Vsin/ = 120 + 30 km/s 
(Boehm & Catala [19951 . 

The DOG-2 wavelet transform (Figure |6j, which is an ex- 
cellent tool for studies of chaotic noise, shows a different pic- 
ture from the Morlet-6 transform: We see individual brightening 
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Fig. 7. Upper panel: The scalegram for the 2 1 days of the MOST 
data for HD 37806 based on the DOG-2 wavelet (continuous 
line) for the data were sampled at the MOST orbital period of 
0.07 d. The horizontal axis gives the time scale in days. Because 
only slope values and their changes are important here, the ver- 
tical scale has an arbitrary zero point. The broken line has slope 
3 in log-log units. The dotted line gives the scalegram for the 
shortest time scales, based on 9.5 days of data sampled at one- 
minute intervals (see Figure [8]); it has been shifted vertically to 
match the results for the coarsely sampled data. The arrow points 
at time scales where visual inspection suggests appearance of co- 
herent structures in the DOG-2 wavelet transform. Lower panel: 
The continuous line gives the median standard deviation <x* for 
the full set of the MOST data (sampled at one minute), calcu- 
lated for progressively longer time scales and corrected for the 
observational random error of <x s = 0.0028, as described in the 
text. The dotted and broken lines give the ESO and the ASAS 
median dispersions at a given scale, corrected for observational 
errors ctr of 0.007 and 0.017, respectively. 



events which become larger for longer time scales. Some events 
die out quickly or may be just observational noise effects, but 
some grow in brightness for longer time scales. The growing 
strength with temporal scale is a characteristic feature for a self- 
similar noise process. This is discussed in the next section. 

3.3. Brightness variability as stochastic noise 

An underlying stochastic process is considered to be self-similar 
with similarity exponent a when it behaves as: X(At) = A a X(t). 
Scargle et al. (1993) considered noise characteristics of accre- 
tion sources and introduced the concept of a "scalegram" func- 
tion to measure the self-similarity properties of such processes. 
The same approach with an extensive discussion and utilization 
of wavelets for cataclysmic variable data appeared in Fritz & 
Bruch ( 1998). Indeed, the scalegram is particularly conveniently 
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implemented using the wavelet transform results: In practice, it 
is the wavelet-summed power for progressively increased time 
scales as they change along the vertical axis in Figs. [5] or [6] 
With N points of input data, for the scale (number of equidis- 
tant data intervals) incremented as S = 2,4, 8, the scalegram 
V(S) for the scale S is denned as: V(S) oc (2S/N) x £ Wf(S), 
where W/{S ) are the individual wavelet transform components 
for a given scale. If V(S ) shows a linear behaviour in logarithmic 
units corresponding to V(S) oc S 2a , then the underlying process 
is considered to be self-similar with exponent a. 

For HD 37806, the DOG-2 wavelet transform appears to pro- 
vide the best description of isolated and/or superimposed flares. 
In Fig. [6] we directly see the intensifying wavelet strength in the 
upward direction, towards longer scales. When expressed as the 
scalegram in log-log units (Figure [7] the upper panel), the pro- 
gression takes the shape of a linear dependence on time scales 
ranging approximately between 0.25 days and 4 days. For the 
MOST-orbit averaged data, the shorter time scales (shorter than 
the MOST orbit length) are not accessible, while scales longer 
than about 4-5 days become affected by edge effects of the 
full data sequence. Within the accessible range, the logarithmic 
slope appears to be about 3 (not a formal fit, just a rough approx- 
imation), so that the noise process is characterized by the expo- 
nent a =i 1.5. The sign and the magnitude of the exponent mean 
that longer noise events are stronger in intensity. Of course, any 
process with a > has such a property, but we note that this 
particular value of the exponent is identical to that seen by Fritz 
& Bruch (1998) in many cataclysmic variables. However, for 
HD 37806, the time scales involved are hours and days rather 
than minutes as observed in CV's. 

The MOST-orbit frequency sampling of the data at 0.07 d 
sets an obvious limitation for the shortest time scales. This could 
be partly remedied by using the second half, 9.5 days long por- 
tion of the MOST observations which were obtained at one 
minute intervals without repeated scattered-light gaps. Still, even 
for this sequence, some data dropouts were present resulting in 
97% temporal coverage; those single, occasional gaps were filled 
by spline interpolated points. The DOG-2 wavelet transform 
looks very different without any obvious indication of growing, 
coherent structures (Figure [8). The small features which we do 
see can be explained by minor instrumental discontinuities and 
other data imperfections which became relatively more impor- 
tant when the fluctuations in stellar flux were very small. 

The wavelet transform and the related scalegram are excel- 
lent tools in defining the temporal scales and the noise-strength 
progression with scale. But the size or strength of the fluctua- 
tions is difficult to gauge. To relate the scalegram results to the 
observed data scatter measured by the standard mean error o~ b S , 
a simple analysis of the scatter was done for the whole available 
MOST dataset: Values of the mean standard error cr ts were cal- 
culated along the time sequence for increased scales (intervals) 
starting at single minutes and ending at 10 days, making sure 
that each interval had a meaningful number of at least 5 points. 
Then, for each scale S, the median cr mei i{S) value was found 
from all individual values of cr ois as a characteristic noise num- 
ber for that scale. Although very simple, this approach has led 
to surprisingly well defined values of the typical fluctuations as 
they varied with the length of the interval. The results are plot- 
ted in the lower panel of Figure [7J after being corrected for the 
observational random error ctr as cr+(S) = o~ 2 med {S) - cr R . The 
assumption for the MOST data was that <x R = 0.0028 as esti- 
mated for the shortest scales and in agreement with the expected 
observational random error at V - 8. The results confirm the 
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Fig. 8. The DOG-2 wavelet applied to the HD 37806 data sam- 
pled at one minute intervals during the second half of the MOST 
observations. Note the qualitatively different picture from that in 
Figure [6] 



scalegram progression and may even confirm that a change in 
the noise characteristics was present around 0.25 days. 

The same procedure of the cr\ determination as for the 
MOST data was repeated for the longer but more poorly sam- 
pled ESO data (Section [23] and Fig. [D and the ASAS data 
(Section [2.41 and Fig. |2}. The MOST and the ground-based re- 
sults appear to agree very well and overlap in Figure [JJ for the 
stated observational error of the ESO observations of <x s = 
0.007. For the ASAS observations the error has been assumed to 
be ctr = 0.017, which is a less well defined number but agrees 
with the typical accuracy of 0.015 to 0.02 mag. The brightness 
fluctuations of HD 37806 appear to grow in size up to the ac- 
cessible scales of some two thousand days, confirming the "red- 
noise" characteristics of the variability. 

4. Discussion and interpretation 

Because of the lack of spectroscopic support, our data cannot 
unequivocally determine the cause of variability of HD 37806. 
Some guidance can be found in extensive recent results for 
T Tauri stars to which Herbig Ae stars appear to be more similar 
in terms of their activity than to Herbig Be stars (Mottram et al. 
(2007). Of various processes identified in T Tauri stars (Bertout 
119891 Bertout et al. 119881 Bouvier et al. 120071 1. one may con- 
sider: (1) photospheric spot rotational modulation, (2) variable 
obscuration by the dusty accretion disk, (3) instabilities in the 
accretion shock funnel on the stellar surface, (4) instabilities in 
the accretion disk. In addition, Herbig Ae stars appear within the 
6 Set instability strip so that pulsations are also possible. 

The wavelet analysis tells us that variability of HD 37806 
takes the form of upward-directed spikes, rather than dimming 
events. This suggests that it is caused by accretion instabilities 
on the surface or in the space adjacent to the star. Rich and com- 
plex accretion phenomena in young stars have been thoroughly 
described in the monograph of Hartmann (2009). While consid- 
ered by us as the most likely process from the start, the accre- 
tion explanation has recently found a very strong support in the 
discovery of Wheelwright et al. (2010) that the star is a binary 
system. This discovery needs confirmation as the angular sep- 
aration is relatively small for the method used, 0.025 + 0.001 
arcsec, and the data for HD 37806 are of moderate quality. The 
implied physical separation of components for the distance of 
300 - 500 pc is of the order of 0.7 - 1.0 AU indicating that the 
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accretion process involves relatively large spatial dimensions; 
this would agree with our results that the red-noise variability 
extends to temporal time scales as long as several years (see be- 
low). 

An explanation of the observed variability of HD 37806 by 
stellar pulsations does not seem to be viable. Light variations of 
HD 37806 in no way resemble stellar pulsations observed for 
another Herbig Ae star HD 142666, which was also observed 
by MOST (Zwintz et al. 12009b . Also, HD 37806 is hotter than 
the hot edge of the classical Delta Scuti instability strip. Still, 
we carefully analyzed the MOST data for HD 37806 in the fre- 
quency range where 5 Scuti-like pulsations would be expected, 
i.e. for frequencies between 5 and 80 c/d or periods between 
about 5 hours and 20 minutes. The only detectable, significant 
peaks in the amplitude spectrum can be attributed to the instru- 
mentally induced frequencies which are well characterized for 
the MOST data. 

The quasi-periodic signal with a period of about 1 .5 d, which 
is visible in the Fourier (Figure|4| as well as in the Morlet trans- 
form (Figure |5j, finds the simplest explanation in an accretion 
area periodically brought into view by rotation of the star. Note, 
that this appears to be a bright, weakly-defined signature rather 
than one corresponding to a dark, sunspot-like feature. The ob- 
served Vsin/ = 120 ± 30 km/s (Boehm & Catala [1995T l can 
be reproduced for this period - assuming rigid rotation - by 
R = (3.5 ± 0.8 ^ )/ sin i. 

The main finding of this work of self-similar brightening 
events combining into chaotic red-noise (Figures [6] and [7]) may 
be interpreted as manifestations of instabilities in an accretion 
disk around the star with progressively larger areas of the accre- 
tion disk producing longer-lasting events. The instabilities pro- 
duce undetectable (<0.0003 of the mean flux) variations for time 
scales of minutes, become very well defined (at about 0.003 of 
the mean flux) for time scales of hours, grow to 0.010 of the 
mean flux at time scales of 1 - 10 days and continue increasing 
for time scales of weeks to a few years where they reach 0.025 
of the mean flux (see Figure |7}. Although these variations are 
very well visible and characterized thanks to the good quality of 
the observations, they remain relatively small, indicating that the 
star dominates in the mean brightness. 

If interpreted by circumstellar events with time scales dy- 
namically related to sizes, for a 3 M G star, dimensions corre- 
sponding to scales of hours to years would range between a 
fraction of a stellar radius to ^ 3 AU. Such a range of dimen- 
sions may correspond to an accretion disk similar to those ob- 
served in T Tauri stars or one created by the matter falling from 
a close companion. The interferometric observations of Eisner et 
al. (|2004) did not rule out the presence of a binary companion at 

< 1 AU and were unable to determine the inclination of the cir- 
cumstellar material around the star. The new spectro-astrometric 
method of Wheelwright et al. (120101) did detect a companion at 

< 1 AU, but the detection is rather preliminary and requires con- 
firmation. 

Does the red-noise variability keep on growing at time scales 
longer than 9 years covered by the ASAS (Section l2~4l i project? 
A comparison of V magnitudes obtained from various sources 
(Section 12.11 and extending over about half a century, is shown 
in Figure [9] Although the data are of a somewhat non-uniform 
quality, there is a clear indication that variations of HD 37806 are 
apparently small, not exceeding +0.05 of the mean flux, which is 
not much more than accuracy of some of the early observations. 
The small spread in the measured V magnitudes is in fact sur- 
prising taking into account that some of the archival observations 
were transformed from the y-band of the ubvy photometry and 
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Fig. 9. V-band magnitudes of HD 37806 from the literature 
(Hardie et al. 119641 Warren & Hesser 119771 de Wi nter e t al. 
I200TI see Section ITTl crosses), ESO (Manfroid et al. [19"9"T] see 
Section 1231 filled circles) and ASAS (Section l2~4l open circles) 
with sizes of symbols logarithmically scaled depending on the 
number of data points in a given season (compare with Figures [TJ 
and|2]i. While the literature plotted data correspond to typically 2 
- 5 observations obtained within a few days, the ESO and ASAS 
observations correspond up to 80 observations per season. The 
large cross marks the mean Tycho-2 V magnitude with its uncer- 
tainty assigned here for the duration of the Hipparcos mission. 



that the star is an emission-line object; both circumstances might 
have easily increased the observational scatter through problems 
of photometric transformations. We note that the only indication 
of systematic changes is the presence of a small brightness de- 
crease which may have taken place sometime between 1995 and 
2000. 

5. Summary 

Although the Herbig Ae/Be star HD 37806 was a secondary ob- 
ject during the MOST observations of the Sigma Ori field, the 
results may be of primary importance to studies of variability of 
these young, massive stars. Variability of HD 37806 appears to 
be similar in general characteristics to that of the T Tauri star 
TW Hya (Rucinski et al. (2008) in that both stars - observed by 
MOST and analyzed in a similar way - showed brightness fluc- 
tuations apparently due to instabilities in the matter accretion 
process. An accretion disk may be fed by circumstellar matter or 

- in the case of HD 37806 - most likely by a companion. The 
process is unstable and produces progressively larger variability 
amplitudes at longer temporal scales, a property which can be 
characterized as "red noise". 

The variations in TW Hya and HD 37806 differ in terms 
of temporal as well as magnitude of variations. For TW Hya, 
the red noise contained discrete, relatively large (typically 0.1 

- 0.2 of the mean flux), periodic variations lasting typically a 
couple of weeks. They evolved in their period-length towards 
shorter periods (a well defined variation evolved from about 6 
days to 3 days in three weeks and another appeared at about 
7 days). For HD 37806, variations are much smaller (typically 
0.01 - 0.02 of the mean flux) which may reflect a much larger 
brightness of the star itself, and we see no indications of pe- 
riodic components similarly drifting in period: A single, weak, 
periodic signature of about 1.5 days did not show any temporal 
evolution in time and may be related to inhomogeneities on this 
rapidly rotating star. However - importantly - brightening events 
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of HD 37806 showed clear self-similar characteristics, growing 
in intensity with time scale. The self-similar characteristics of 
stochastic variability were established for the time scale range of 
about 4 hours to 4 days. The self-similarity exponent (a 1.5) 
of the noise was found to be very similar to that observed in cat- 
aclysmic variables at much shorter time scales of a minute. The 
similarity to the CV's may be not accidental in view of the recent 
result of Wheelwright et al. (120101 1 that HD 37806 is most likely 
a binary with the component separation of < 1 AU. 

The archival data on the variability of HD 37806 suggest its 
continuation - with a possible moderation - for long time scales. 
The red noise characteristics do appear in a simple analysis of 
the typical (median) standard dispersion of the noise for progres- 
sively increased time scales. The large temporal and time-range 
extent of the combined MOST, ESO and AS AS data permits one 
to see a well defined progression in the amplitude - scale depen- 
dence up to intervals of a few years. The red-noise fluctuations 
grow for time scales ranging from minutes to a few years from 
below 0.0003 to 0.025 of the mean flux. The fragmentary and - 
by their nature - somewhat non-uniform archival V magnitude 
data do not suggest any drastic noise increase (< 0.05) for still 
longer time scales extending to half of a century. 
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